Suppression of superconductivity at nematic critical point in underdoped cuprates 
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A nematic quantum critical point is anticipated to exist in the superconducting dome of some 
high-temperature superconductors. The nematic order competes with the superconducting order 
and hence reduces the superconducting condensate at T = 0. Moreover, the critical fluctuations of 
nematic order can excite more nodal quasiparticles out of the condensate. We address these two 
effects within an effective field theory and show that superfluid density p"{T) and superconducting 
temperature Tc are both suppressed strongly by the critical fluctuations. The strong suppression of 
superconductivity provides a possible way to determine the nematic quantum critical point. 
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I. INTRODUCTION 

The strong electron correlation in high-temperature 
superconductors (HTSC) is able to drive an electronic 
nematic phaseA^-i, which preserves translational symme- 
try but breaks rotational symmetry. In the past decade, 
there have been a number of experimental signatures 
pointing to the presence of nematic ordering transition 
in some HTSCs^"— . On the basis of these experiments, a 
zero-temperature nematic quantum critical point (QCP) 
is supposed to exist at certain doping concentration Xc 
in the superconducting (SC) dome^^^^—. Generally, the 
nematic order has two impacts on the SC state. First, 
it competes with the SC order—"—. Second, the nematic 
order parameter couples to the gapless nodal quasiparti- 
cles (QPs), which are believed to be the most important 
fermionic excitations in unconventional superconductors 
with (1^2 _y2 energy gap. The latter coupling is singular 
at the nematic QCP Xc, and has stimulated considerable 
theoretical efforts^"—. A recent renormalization group 
analysisii showed that it leads to a novel fixed point at 
which the ratio between gap velocity wa and Fermi ve- 
locity vp of nodal QPs flows to zero, v\/vf — > 0. 

Although a zero-temperature nematic QCP is expected 
to exist somewhere in the SC dome^i^rJ^, shown schemat- 
ically in Fig. ([T]), its precise position, and even its very 
existence, has not been unambiguously confirmed by ex- 
periments so far. It is therefore always interesting to seek 
evidence which can help convincingly confirm or discon- 
firm the existence of such point. In this paper, we study 
the superfiuid density p'^(T) and the SC temperature Tc 
at the supposed nematic QCP Xc- If P^{T) and Tc exhibit 
sharply distinct behaviors at Xc, then the nematic QCP 
may be detected by measuring these quantities. 

HTSCs are known to be doped Mott insulators, so their 
superfluid density is much smaller than that of conven- 
tional metal superconductors. At T = 0, the superfluid 
density in underdoping region depend o^^i^^ linearly on 
doping X as p*(0) — x/a?, where a is the lattice spac- 
ing. At flnite T, certain amount of nodal QPs are ther- 
mally excited out of the SC condensate. Lee and Wen ar- 
gued that these normal nodal QPs can efficiently deplete 



the superfluid densityi^. Formally, the superfluid density 
contains two terms, p'{T) = p"(0) - /3"(r), where p"(r) 
is the normal QPs density. Setting p^{Tc) — allows for 
an estimate of the critical temperature T^. Employing 
a phenomenological approach, Lee and Wen^* obtained 
Tc oc /O*(0) cx reproducing the Uemura plo1>i^. 

Once a nematic ordering transition occurs at Xc, the 
superfluid density and Tc will be substantially changed. 
As va/vf due to the critical nematic fluctuations, 
it seems that Tc — >■ 0, i.e., superconductivity would be 
completely suppressed at Xc- This argument is certainly 
oversimplified since the above expression of Tc is obtained 
in the non-interacting limit. However, this qualitative 
analysis does indicate the importance of the critical ne- 
matic fluctuations, and indeed motivates us to perform a 
quantitative computation of the renormalized p'^ (T) and 
Tc after taking into account the nematic fluctuations. 

The nematic order affects p^{T) in two ways. On the 
one hand, since the nematic order competes with the 
SC order, it reduces ^"(0). This reduction can be ex- 
amined by studying the competitive interaction between 
nematic and SC order parameters. On the other, the 
critical nematic fluctuations can excite more nodal QPs 
out of the condensate, compared with the case without 
nematic order. As a consequence, p"(T) is enhanced and 
the superfluid density is further suppressed. We shall ac- 
cess this effect by generalizing the phenomenological ap- 
proach proposed in Refii^. The velocity anisotropy plays 
an essential role in these considerations. After explicit 
calculations, we find that superfluid density p'^iT) and 
Tc are both significantly reduced due to critical nematic 
fluctuations, indicating a strong suppression of supercon- 
ductivity at nematic QCP Xc (see Fig. ([T|)). 

The rest of the paper is organized as follows. In Sec. II, 
we address the competitive interaction between SC and 
nematic order parameters and calculate zero-T superfluid 
density. In Sec. Ill, we calculate the density of normal 
QPs after taking into account fermion velocity renor- 
malization due to critical nematic fluctuation. Based on 
these calculations, we predict a strong suppression of su- 
perconductivity at nematic QCP. In Sec. IV, we present 
a brief summary of our results, and also discuss the pos- 
sible experimental determination of the nematic QCP. 



2 



A 

\ \ 






\ \ 
\ \ 
\ \ 

\ \ 

\ \ 

\Nematic\ 

\ \ 






\ \ 

afm\ / 


1 

. sc 


— > 



X 



FIG. 1: Schematic phase diagram of HTSC. The suppression 
of Tc at nematic QCP Xc is represented by the dashed Une. 



II. COMPETING ORDERS AND 
ZERO-TEMPERATURE SUPERFLUID DENSITY 

We first consider the renormalized zero-T superfluid 
density at nematic QCP. Based on phenomenological 
grounds, we write down a free energy of two competing 
orders, 



F = 



1 
2m 



F^ + F 



i'<f> 



F,. 



(1) 



where ^ and (p are the SC and nematic order parameters, 
respectively. The strength of the competitive interaction 
between ip and (j) is represented by a positive constant 7. 
Such type of effective model has been adopted to describe 
competing orders in various superconductors^SiSi. 

In the absence of nematic order, the mean value of tp 
is lipl = ^/aj]3. To be consistent with experiments, the 
parameters a and /? must be properly chosen such that 
4|?/)p = 4q;//3 = p*(0) — xja}. In the presence of nematic 
order, \\p\ will be renormalized by the ^i\P'(p?' term. The 
quantum fluctuation of nematic order (p is very strong 
and actually singular at nematic QCP Xc, so <p should 
be regarded as a quantum-mechanical field. However, we 
can consider -0 as a classical field and ignore its quantum 
fluctuations, provided Xc is not close to the SC QCP xq. 

The free energy of nematic order, i^^, should be spec- 
ified now. Analogous to the free energy of SC order ip, 
F^ contains a quadratic term (p^ and a quartic term cp'^ . 
However, the additional coupling between nematic order 
and nodal QPs introduces an extra term. The action 
describing this coupling is given by^rii 



S-i!^ — 



(2^ 



-f^'^^(-iw + VFk2T'- + VAklT'')'i'2i], (2) 
J d^xdT[XocP{^lT^^u + ■^lT^'^^2^)], (3) 



where r^'^'^ are Pauli matrices and the flavor index i 
sums up 1 to N. 5'i represents nodal QPs excited from 
(f , f ) and (-f , -f ) points, and *2 the other two. The 
effective action of 6 has the fom>ii 



Sd, 



d^ 



{2n) 



l{q' + r + n(g)) c^' + 



(4) 



which is obtained by integrating out A''- flavor nodal QPs. 
Here, r is the tuning parameter for the nematic ordering 
transition, with r — defining the nematic QCP Xc- The 
polarization function Tl{q) comes from nodal QPs, and, 
to the leading order of 1 /TV-expansion, is defined as 



n(e,q) 



where 



N 



(27r)3 



(5) 



Go(c^,k) 



(6) 



is the free propagator for nodal QPs (free propagator 
for nodal QPs can be similarly written down). The 
polarization function H(e, q) has already been calculated 



previousl; 



n(e,q) 



,11.15 



and is known to have the for 



N 



(62 



(7) 



Note there is no direct interaction between ip and nodal 
QPs ^. Indeed, are excited on top of a SC order ip. 
Moreover, they have sharp peak and long lifetime in the 
SC dome in the absence of competing ordersii, so their 
coupling to Ip must be quite weak. 

Starting from Eq. ([T|) and Eq. (jj]), we can compute 
the correction to zero-T superfluid density due to the 
competition between SC and nematic orders. To this 
end, we need to minimize the effective potential ofSC 
order, V[ip]. The bare potential for ip is simply Vo[ip] = 
-aip^ + ^ip^. It receives an additional term Vi [ip] due to 
the SC-nematic competition. This additional term will 
be calculated using the methods presented in Refs. — 
The corresponding partition function is 



zmr)] 



J 'Dip{r, t) exp 



Sd, 



(8) 



where — J d'^rF^. The saddle-point equation for ip 
reads 



5lnZ[iJj{r)] 
6ip{r) 



0, 



which gives rise to 



-a + /3V^2(r)--V^+7/[V'] 

TO 



(9) 



■p{r) = 0. (10) 
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Here, is the expectation value of ({P' . At the one-loop 
level, it has the form 



1 



(27r)3 



U{q) + 



2 ' 



(11) 



where Il{q) and 7'!/'^ represent the contributions due to 
nodal QPs and SC order, respectively. Although the 
q^ term appearing in the denominator of f[ip] is much 
smaller than Il{q) at low energy, it cannot be simply ne- 
glected since f[ip] would be divergent without such term. 
It is now straightforward to get 



svM 

Sip 



= 27^/ 



which then leads to a renormalized potential 



(12) 



(13) 



To calculate the renormalized zero-T superfluid density, 
Pp(0), one needs to minimize the effective potential V[tp] 
by taking 



syjip] 
5ij 



0. 



(14) 



The renormalized and therefore /5f((0) can be ob- 
tained from the solution of the following equation 



aip + /3V^ -H7'(/'/[V'] = 0. 



(15) 



In order to see the role of nodal QPs, we first ignore 
the fermion contribution in f[ip] by assuming Tl{q) — 0. 
In this case, the integration over q in /[ip] can be exactly 
performed, yielding 



is expected and trivial. As 7 increases continuously, 
and p^{0) both decrease rapidly. To estimate this effect 
more quantitatively, we assume an ultraviolet cutoff A = 
lOeV, and choose a = 2.5 x 10~3eV. Moreover, we con- 
sider a representative bare velocity ratio v\/vf ~ 0.075, 
which is an appropriate value for YBa2Cu306+#^. We 
further choose p^{0)/m = Aa/m/3 = 10~^eV, which cor- 
responds to Tc w 20K. It is now useful to introduce a 
parameter given by 70 = P/^oCe with = l/2ma, and 
define a dimensionless coupling constant 7/70- Fig. ^a.) 
shows that pf^ (0) is strongly suppressed by the nematic 
order and is completely destroyed when 7 is large enough. 
It is also obvious that bare velocity ratio plays a vital 
role: a larger anisotropy causes smaller drop of p''(0). 
This property provides further evidence that nodal QPs 
can not be simply ignored. 




FIG. 2: (a) pji{0) / p" (0) for four different values of bare ratio 
va/vf; (b) Tc/Tc- The open circles etc. are guides to eyes. 



= — (A-V7^). 
The renormalized potential becomes 

zTT ovr 



(16) 



(17) 



The cubic term induced by nematic order turns the SC 
transition to first orde r^^'^'^ , with critical point being 
ac = ^ — Q^^ip ■ However, since the gapless nodal QPs 
are present even at the lowest energy, the polarization 
function H(g) should be included in the effective action 
S(t). After including the polarization Tl{q) into f[ijj], in- 
tegration over q can not be carried out analytically, and 
numerical method will be used. After numerically solving 
Eq. (|15|) . we found a critical value 7c for the competi- 
tive interaction between SC and nematic orders. The SC 
transition remains continuous if 7 < 7c, and is driven to 
first order when 7 > 7c. 

We now discuss the effects of the competing nematic 
order on the zero-T superfluid density. When 7 is zero or 
very small, the nematic order does not change \ip\'^ , which 



III. VELOCITY RENORMALIZATION AND 
SUPPRESSION OF SUPERCONDUCTIVITY 

We then turn to calculate the density of normal nodal 
QPs, p"(T). We first briefly outline the phenomenologi- 
cal approach of Lee and Wen^^ and then generalize it to 
the case with critical nematic fluctuation. In the scenario 
of Lee and Wen, a crucial problem is how to assess the 
roles played by nodal and antinodal QPs in the destruc- 
tion of SC condensate. As revealed clearly by extensive 
experiments^^, a pseudogap phase exists above Tc. The 
pseudogap turns out to have the same d-wave symme- 
try as the SC gap^^, which implies the antinodal QPs 
remain gapped as temperature increases across Tc. It 
was argued^^ that the SC state is primarily destroyed by 
the thermal proliferation of the low lying, gapless nodal 
QPs, and the antinodal QPs are only spectators. Lee and 
Wen further assume that the nodal QPs do not carry 
superflowi^, so the fermion spectrum is shifted by the 
vector potential to -E(k, A) = -E(k) -I- -Vfe • A, where the 
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QPs energy is linearized as E(k) = \J v\k\ + v\k\ near 
nodes and vj, is the normal state velocity. The 

electric current is — — p^^^A^, where the superfluid 

and the 

normal QPs density is derived from the free energy^ 



tensor can be written as p^y(r) — p*(0)(5, 
The zero-T superfluid density is /9*(0) — x/c 



F(A,r) = -T^ln(l 



-£;(k,A)/T 



(18) 



by the formula j^p^AT) = -'^Y.^^r^TrM- the 
non-interacting case, the fermion velocities are constants, 
so one gets 



1 c ,„s X 2 In 2 Up „ 
m ma tt va 



(19) 



which exhibits a linear temperature dependence, in agree- 
ment with experiments^. At T = 0, p''{T) takes its max- 
imum value, /3'*(0). As T is increasing from zero, p^{T) 
decreases rapidly due to thermally excited nodal QPs and 
eventually vanishes as T Tc- By taking p^{Tc) = 0, it 
is easy to obtain 



Tc oc /9^(0) cx 



vp ma? 



(20) 



This linear doping dependence of is well consistent 
with the Uemura plot. The same results were later re- 
produced by means of Green's function method^^. The 
Fermi liquid (FL) corrections to these results were also 
investigated22i2^. An important fact is that both p'^{T) 
and Tc contain the velocity ratio v/^/vp, so these two 
quantities are expected to be significantly affected by the 
critical nematic fluctuation which can cause a nontrivial 
renormalization of the velocity ratio. 

The above approach of computing Tc is applicable 
when the nodal QPs are well defined. At the nematic 
QCP Xc, the critical fluctuation of nematic order couples 
strongly to gapless nodal QPs and may lead to break- 
down of FL behavior. As a consequence, the nodal QPs 
might no longer be well defined QPs. In a strict sense, 
the validity of the simple d-wave BCS theory^ and its 
FL-interaction generalizations22i2^ are both doubtful. In 
order to proceed, here we assume that the basic approach 
of Refiii^ is still valid after the free energy of QPs receives 
a singular correction due to nematic fluctuation. 

The free fermion propagator shown in Eq. (|6]) is renor- 
malized by nematic order to 



G(w,k) 



1 



Go\uj,k) - S(w,k)^ 
where the self-energy is 
detPq 



S(c^,k) 



{2ny 



-Go(a; + e,k + q) 



(21) 



(22) 



The pole in G(w, k) determines the renormalized energy 
E{uj, k) = E{uj, k) -I- 5E{lo, k), where E{uj, k) is given by 



pole of free propagator Go{uj, k) and SE{ijj, k) represents 
the interaction correction to energy. To describe this 
renormalization procedure, it is convenient to expand the 
self-energy formally as 

T.{uj, k) = -iwEo + Siuj^fcir^ -I- ^2VAk2T'' , (23) 

which exhibits unusual logarithmic behaviors since 
^i,2{k) contains such a term as ln(A/fc). We combine 
this self-energy with G(^^(fc), and find that the bare ve- 
locities acquire strong fc-dependence, vf,a — >■ vp^A{k)- 
VF,A{k) are very complicated functions of fc, and thus 
can not be written down analytically. After straightfor- 
ward numerical computation, we show the fc-dependence 
of fermion velocities and their ratio in Fig. ([3|), for the 
representative bare value va/vf — 0.075. It is clear 
that both vpik) and VA{k) vanish as /c — >■ 0. Neverthe- 
less, vpik) approaches zero much more slowly than VAik) 
does. Therefore, the velocity ratio vanishes at the lowest 
energy, va/vf 0, which recovers the extreme veloc- 
ity anisotropyii. Such velocity renormalization causes 
breakdown of FL behavior^^'^^ and can faithfully char- 
acterize the non-FL nature of QPs at Xc- 
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FIG. 3; Flow of uf, i^a, and va/vf with momentum k for 
bare ratio va/vf = 0.075. As fc — > 0, both velocities are 
reduced down to zero. However, va decreases much faster 
than vf, so velocity ratio va/vf 0. 

Although the nodal QPs are strongly damped by the 
nematic fluctuations and thus no longer well-defined at 
Xc, it is still reasonable to assume that they do not carry 
superflowiS, which allows us to write -E(k, A) = E{\i) + 
-Vfe • A, where the normal velocity becomes strongly 
fc-dependent. Now the new free energy is 



F(A,r) 



-T^ln(l 

k.cr 



,-£;(k,A)/TN 



(24) 



from which we obtain a renormalized superfluid density 
pkiT) = Pr(0) - pUn (25) 

rf2k 4(^yy.|(fc)fe?+.i(fc)fci/T 



4 
T 



-^26) 
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At nematic QCP, Tc is changed from its original value to 
a renormalized value T'^. When T — T^, the renormalized 
superfluid density vanishes, p^{T^) — 0, so that 



(27) 



which builds a relationship between and Tc . The ratios 
of p^{T) / p" (T) and T'^/T^ can be numerically calculated 
according to the above three equations. To include the 
influence of interaction corrections, the velocities vf,a 
appearing in pf^(O) are also replaced by VF,A{k). The 
numerical results of T'^/Tc are shown in Fig. [2lb), and 
those of p^{T)/ p'^{T) given in Fig. (g]). The suppression 
of superconductivity is prominent even when the compet- 
itive interaction between SC and nematic orders is quite 
weak, and is further enhanced by this interaction. As 
aforementioned, the SC phase transition at remains 
continuous as long as 7 is small. 
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FIG. 4: Superfluid density at finite temperature, p^(T). (a): 
7 = 0; (b): 7/70 = 10"^ 

Clearly, the extent to which p'^{T) and Tc are sup- 
pressed by the nematic order depends on a number of 
parameters, such as a, 7, and va/vf- It is helpful to 
make more quantitative analysis. The magnitudes of a 
as well as A have already been chosen. The bare velocity 
ratio is still fixed at va/vf — 0.075. Fig. [2][b) shows 
that Tc is roughly suppressed by 25% when 7 = 0, and 
by 50% when 7/70 « 1.5 x lO^^. The drop of p'{T) 
can be analyzed in a similar way. An important impli- 
cation of Fig. (HI is that p^ (T) is more strongly reduced 
at higher temperatures. Experimentally, the suppression 
of Tc can be tested by measuring resistivity and Meiss- 
ner effect, and the suppression of superfluid density may 
be probed by measuring London penetration length Al, 
preferably using microwave techniques2^i2£ based on the 
relationship p%T) cx A£^(r). 

The strong suppression effects rely on the divergence 
of correlation length ^ of the nematic field </> at the QCP 
Xc- As one moves away from the point Xc, the strong 
suppression of superfluid density and Tc will be rapidly 
diminished. In fact, once correlation length ^ becomes 



flnite, the interaction between nematic fluctuation and 
nodal QPs is no longer singular. Consequently, the self- 
energy of nodal QPs, does not display logarith- 
mic behavior, which indicates that the fermion velocities 
only receive unimportant renormalizations. In this case, 
the density of normal QPs is not significantly modified 
by the competing nematic order, since the velocity ratio 
does not depart far from its bare value. The influence 
of nemtic order on /9^(0) is also weakened. For x > Xc, 
the fluctuation of </> is rather weak and does not change 
p'^{0) much if 7 is not very large. For xq < x < Xc, the 
fluctuation of (j) is gapped and also leads to much smaller 
change of (0) compared with what happens at nematic 
QCP Xc- Therefore, the Tc curve generally follows the 
Uemura plot in the whole underdoped region, but devi- 
ates from the linear behavior in the close vicinity of Xc- 
Since the suppression of superconductivity is most pro- 
nounced at Xc, it may help to flnd the position of the 
predicted nematic QCP if such a point really exists. 



IV. SUMMARY AND DISCUSSION 

In summary, we have considered the effects of nematic 
fluctuation on superfluid density and critical temperature 
Tc in the vicinity of nematic QCP in the contexts of some 
HTSCs. On one hand, the critical fluctuation of nematic 
order parameter reduces zero-T superfluid density as a 
result of ordering competition. On the other hand, it cou- 
ples strongly to the gapless nodal QPs and excites more 
normal QPs out of the SC condensate by triggering an 
extreme fermion velocity anisotropy. Both of these two 
effects combine to signiflcantly suppress superconductiv- 
ity. Therefore, we have predicted a dip shape reduction 
of Tc at nematic QCP Xc, which is schematically shown 
in Fig. (P. 

The dip shape of Tc reminds us of 1 /8 anomaly^r— 1^"— , 
which refers to an anomalous suppression of supercon- 
ductivity at doping a; = 0.125. A sudden drop of 
p^(0) was also observed^^ in La2-a;Sra;Cu04 (LSCO) at 
X = 0.125. It appears that p'^{T) and Tc near the ne- 
matic QCP bear some similarity to the basic features 
of these experiments. However, we refrain from identi- 
fying the 1/8 doping as the anticipated nematic QCP 
for several reasons: i) 1/8 anomaly is usually attributed 
to the formation of static stripe order-Si^Si; ii) it is not 
clear why nematic QCP exists precisely at x = 0.125, 
not elsewhere; iii) low-T dc thermal conductivity k/T 
was predicted to be enhanced at nematic QCP^"^, whereas 
an early transport measurement found a drop of k/T in 
LSCO at a; = 0.125^. 

Recently the critical nematic fluctuations are shown 
to induce unusual behaviors in several quantities, 
such as QPs spectral functioiii^, specific healii^, nu- 
clear relaxation rate^^, thermal conductivity^, and QP 
interferencei^. Normally, a single observation alone is 
not able to uniquely fix the nematic QCP. Should all or 
most of these predictions, including the suppression of 
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superconductivity, be observed, the nematic QCP might 
be determined. However, the absence of a sharp nematic 
QCP does not necessarily mean the absence of nematic 
phase, since the sharp QCP may be rounded by disorders 
and become a crossover^°. In case this smearing occurs, 
many of the anomalous critical behaviors are weakened, 
or even destroyed, and the nematic phase most possibly 
shows its existence in the uj <^ T regime. 

In the calculations presented in this paper, we have as- 
sumed that the SC order parameter -0 is a classical field 
and neglected its quantum fluctuation. This assumption 
is expected to be a good one if the nematic QCP Xc is 
deep inside the SC dome, namely Xc is not close to the SC 



critical point xq- When Xc is close to Xq, the quantum 
fluctuation of SC order may become very important 
and a fully quantum- mechanical treatment will be neces- 
sary. 
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